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PSEUDO-OFFSET MIGRATION: 

BACKGROUND OF THE INVENTION: 

This invention relates to the field of seismic data acquisition and processing. 

Converted-wave (C-wave) processing is more difficult than pure-mode 
processing. Those of the ordinary skill in the art will recognize that converted waves 
result from compressional waves (P-waves), which come into contact with a reflecting 
surface and incite a shear wave (S-wave). The opposite is also true. That is, an S-wave 
may also incite a P-wave. Because P-waves and S-waves have different characteristics, 
converted wave processing, which requires a set of characteristics for the down-going 
wave and a different set of characteristics for the up-going wave, is much more 
complicated than processing only a single set of characteristics for only P-wave data or S- 
wave data. 

For example, accurate migration is critically dependent on both P-wave and S- 
wave velocities. Furthermore, conventional prestack time migration techniques, such as 
the industry standard of Kirchhoff migration, are normally based on the double-square 
root travel time equation. These methods assume a straight ray from the source to the 
scatter point and from the scatter point to the receiver. Other current Kirchhoff prestack 
time migration techniques are often implemented with a travel time approximation (Li et 
al., 2001). The straight ray assumption and the travel time approximation thus, limit the 
ability to produce an accurate subsurface image without a correction for ray bending and 
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travel time errors. For multi-layered v(z) media or vertical transversely isotropic (VTI) 
media, more sophisticated travel time equations are needed. 

Even further, another current prestack migration is equivalent-offset migration. 
Equivalent-offset migration (EOM) is a prestack time migration method proposed by 
5 Bancroft et al. (1998). Equivalent offset migration claims to be computationally efficient 
and weakly velocity-dependent. However, in C-wave processing, equivalent offset 
migration (Wang et al, 1996) depends on the initial velocity even for a single constant- 
velocity flat layer. Thus, EOM is poor for C-wave velocity analysis. 

Finally, many current ocean bottom cable (OBC) techniques require wave- 
io equation datumming. Wave equation datumming adds additional processing steps, 
increases cost, and causes irregular geometry effects on 3-D OBC data processing. 

Thus, there is a long felt need for migration techniques which address one or any 
of the above issues. 
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SUMMARY OF THE INVENTION: 

In one example embodiment of the present invention, a method is provided for 
migrating an input seismic data point having an input source location, an input receiver 
location, and an image location associated therewith. The method comprises determining 
a pseudo-offset, and mapping the seismic data point to the image location based at least 
in part on the pseudo-offset. 

In a further example embodiment of the present invention, a method is provided 
for migration of a seismic data point having an input source location, an input receiver 
location, and a scatter point associated therewith. The method comprises determining a 
projected source location and determining a projected receiver location. The method 
further comprises mapping the seismic data point from the input travel time to a projected 
travel time. The method also comprises determining a pseudo-offset based at least in part 
on the projected travel time, and mapping the data point to the pseudo-offset. 

In an even further example embodiment of the present invention, a method of 
performing velocity analysis on a seismic gather is provided. The method comprises 
computing a moveout travel time based on an initial model and computing an 
intermediate travel time. The method further comprises mapping the gather from the 
moveout travel time to the intermediate travel time, and scanning for velocities based on 
the intermediate travel time. 

In an even further embodiment, a system is provided for migrating an input 
seismic data point having an input source location and an input receiver location, a scatter 
point, and an image location associated therewith. The system comprises means for 
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determining a pseudo-offset, and means for mapping the seismic data point to the image 
location based at least in part on the pseudo-offset. 

In an even further embodiment, a system is provided for migration of a seismic 
data point having an input source location, an input receiver location, and a scatter point 
associated therewith. The system comprises means for determining a projected source 
location, means for determining a projected receiver location, means for mapping the 
seismic data point from an input travel time to a projected travel time, means for 
determining a pseudo-offset based on the projected travel time, and means for mapping 
the data point to the pseudo-offset. 

In an even further embodiment, a system for performing velocity analysis on a 
seismic gather is provided. The system comprises means for computing a moveout travel 
time based on an initial model, means for computing an intermediate travel time, means 
for mapping the gather from the moveout travel time to the intermediate travel time, and 
means for scanning for velocities based on the intermediate travel time. 

In an even further embodiment, a seismic data point is provided. The seismic data 
point is produced by the process of migration comprising determining a pseudo-offset, 
and mapping the seismic data point to the image location based at least in part on the 
pseudo-offset. 

In an even further embodiment, a seismic data point is provided. The seismic data 
point is produced by a method of migration comprising determining a projected source 
location, determining a projected receiver location, mapping the seismic data point from 
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the input travel time to a projected travel time, determining a pseudo-offset based on the 
projected travel time, and mapping the data point to the pseudo-offset. 
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BRIEF DESCRIPTION OF THE FIGURES: 

FIGURE 1 is a diagram outlining an example embodiment of the present 
invention. 

FIGURE 2 is a schematic diagram demonstrating features of an example 
5 embodiment of the present invention. 

FIGURE 3 is a diagram outlining an example embodiment of the present 
invention. 

FIGURE 4 is a schematic diagram demonstrating features of an example 
embodiment of the present invention. 

0 FIGURE 5 is a diagram outlining an example embodiment of the present 

invention. 
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DETAILED DESCRIPTION OF EXAMPLE EMBODIMENTS OF THE 
INVENTION: 

Those of ordinary skill in the art will appreciate that for the purpose of teaching 
the broad aspect of the present invention, example embodiments of the invention are 
described in detail. 

As illustrated by reading Figures 1 and 2 in conjunction, in one example 
embodiment of the present invention, a method of migrating seismic data is provided. In 
alternate embodiments the content and acquisition situation of the seismic data to be 
migrated is land data, dual sensor data, streamer data, VSP data, vertical cable data, or 
any other data from any other acquisition situation that will occur to those of ordinary 
skill in the art. In a further embodiment, the seismic data comprises an input seismic data 
point having an input source location S and an input receiver location R. Those of 
ordinary skill in the art will recognize that such information is often available on the trace 
headers or geometry database. The method comprises determining (1 10) a pseudo-offset 
x and mapping (120) the seismic data point to the image location I based at least in part 
on the pseudo-offset x. In an even further embodiment, the method further comprises 
determining a pseudo-source S* and a pseudo-receiver R*. In still a further embodiment, 
the pseudo-source S* and the pseudo-receiver R* are not collocated. In an even further 
embodiment, the pseudo-offset x comprises a distance between the pseudo-source S* and 
the image location I plus a distance between the pseudo-receiver R* and the image 
location I. In still an even further embodiment, the method further comprises 
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determining a pseudo ray parameter p. In an even further embodiment, determining a 
pseudo ray parameter p further comprises determining a pseudo ray parameter p which is 
constant from the pseudo-source S* to the scatter point D and from the scatter point D to 
the pseudo-receiver R*. In still a further embodiment, determining a pseudo-offset x 
5 further comprises determining a travel time associated with the input source location S 
and the input receiver location R and the scatter point D. In an even further embodiment, 
determining a pseudo-offset x further comprises determining a pseudo travel time from 
the pseudo-source S* to the pseudo-receiver R* via the scatter point D. In an even 
further embodiment, the pseudo travel time and the actual travel time from the input 
10 source location S to the input receiver location R are essentially equivalent. 

In a further example embodiment of the present invention, a method of migrating 
seismic data is developed for multi-layered v(z) isotropic media. To demonstrate this we 
begin with the basic Kirchhoff integral for C-wave prestack time migration. The 
Kirchhoff integral is formulated as the following summation, 

15 output(tQ§)= ^amplitude[filter* input (t)], (1) 

aperture 

where t is the input travel time which is the sum of the P-wave travel time t? from a 
source S to a scatter point D and the S-wave travel time / s from the scatter point D to a 
receiver R as illustrated in Figure 2. The travel time is, 

/ = f P +*s- (2) 
20 t C o is the C-wave vertical travel time which is the sum of the P-wave vertical travel time 
£po and the S-wave vertical travel time * so , 
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'co='po + 'so- (3) 
In one embodiment, for the v(z) media, t ? is computed from the following equation 
(Slotnick 1959), 

«»0*)- f . (4) 

o V 1- ^P v p( r p) 

5 where v P (r P ) is the P-wave interval velocity, and p P is a ray parameter from the source S 
to the scatter point D, which is solved from the following equation, 
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J / 22 =/? P' ( 5 ) 
0 yjl-p?V ? (T ? ) 



where A P is the horizontal distance between the source S and the scatter point D. 
Similarly, in a further embodiment, t s is computed from the following equation, 

ho > 

where v s (r s ) is the S-wave interval velocity, and p s is a ray parameter from the scatter 
point D to the receiver R, which is solved from the following equation , 



JT = TF = s ' (7) 

0 ^~P$Vs( T s) 



where A s is the horizontal distance between the scatter point D and the receiver R. To 
15 replace t? 0 and t$ 0 with fco in equations above, let 

T=T ? +T$, (8) 

and define the interval vertical travel time ratio as 
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ro(r) = dT S /dT ? . (9) 
For the isotropic medium, ^(z) is equal to the interval velocity ratio, that is 

vp(rp)/v s (r s ) = / 0 (r). (10) 

In another example embodiment, we then use the C-wave interval velocity vc(t) defined 
5 byThomsen(1999), 

v g( r ) = v P( r P) rfr P + v s( r s)^s 
dr 

to replace velocities vp(Y P ) and v s (zs). By substituting fo>> *so> v P , and v s with £ C o> vc, and 
yb in the travel time equations (4) - (7), t? and ts is computed using the following 
equations, 

1 
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» sfe ,= f , (14) 

0 f^' Mw 



f: yo<> as) 

0 1 i_n|v^(r) 
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In an even further embodiment, for the isotropic media, we only need to know f C o, 
v c , and to compute f P and f s . Unlike conventional techniques, no effective velocity 
ratio y eff (Thomsen, 1999) is required. In alternate embodiments, the interval vertical 
travel time ratio # is determined directly from the travel time ratio of corresponding C- 
wave and P-wave events such as on stacks or extrapolated from oblique times on prestack 
gathers. 

In an even further embodiment, instead of performing the Kirchhoff summation 
directly as is done in conventional migration, we divide the Kirchhoff summation into 
two steps. The first step is to map each input sample to a new geometry, such that, two 
conditions are met. The first condition is that the travel time from the new source S* 
(called the pseudo-source) to the new receiver R* (called the pseudo-receiver) is 
essentially the same as the travel time computed from the input geometry of the input 
source location S and the input receiver location R, as illustrated in Figure 2. The second 
condition is the ray parameter from the new pseudo-source S* to the scatter point D is 
essentially the same as the ray parameter from the scatter point D to the new pseudo- 
receiver R.* We call the new ray parameter the pseudo ray parameter /?. In a further 
embodiment, the pseudo ray parameter p is computed from the following equations, 

t = t?(P?) + ts<J>s)> (16) 
**(p) = t, (17) 

'*Gp) = 'p(/0 + 's(/>), (18) 
where t* is the new travel time from the pseudo-source to the pseudo-receiver and is 
called the pseudo travel time. As illustrated in Figure 2, the distance between the pseudo- 
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source S* and the image location I plus the distance between the pseudo-receiver R* and 
the image location I is called the pseudo-offset x. In an even further embodiment, with 
the estimated pseudo ray parameter p, the pseudo-offset x is computed as follows 



5 For a given scatter point D, the input samples inside the migration aperture are 
accumulated to different pseudo-offsets x based on the equations (16) - (19). Thus, a 
common-scatter point gather or POM gather is created. We call the first step migration 
mapping. The second step of this migration scheme is to perform moveout correction 
based on equation (18) and stack the POM gather along the pseudo travel time. We call 

10 the second step migration stacking. 

In an even further embodiment, the pseudo-offset migration method is employed 
in the ocean-bottom cable (OBC) acquisition situation. The OBC option assumes the 
elevation difference between shots and receivers and performs the migration without 
wave-equation datuming (Berryhill, 1979). For the OBC data migration, we assume the 

15 velocity function is referenced at the sea level where shots are normally located, so t ? is 
still computed based on equations (12) - (13). However, in an even further embodiment 
of the present invention, the computation for t s is redefined. Assuming d r is the receiver 
depth, v w is the P-wave replacement velocity or the water velocity, and y w is a 
replacement value for the shallow ^ equations (14) and (15) are replaced by 
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d t 




if 



l + ro(^) 




(20) 





P s v w d r 



/>s v cO) 




(21) 



Those of ordinary skill in the art will recognize that in alternate embodiments, 
pseudo-offset migration is still performed based on equations (16) - (18). In an even 
5 further embodiment, the pseudo-offset migration method is employed for OBC, vertical 
cable, VSP, or land data migration or any other acquisition situation that will occur to 
those of ordinary skill in the art. 

Turning now to Figures 3 and 4 in conjunction, a further example embodiment of 
the present invention is presented. A method of migrating a seismic data point having an 

10 input source location, an input receiver location, and a scatter point associated therewith 
is provided. The method comprises determining (310) a projected source location S\ 
determining (320) a projected receiver location R', mapping (330) the input seismic data 
point from the input travel time t to a projected travel time t\ determining (340) a 
pseudo-offset x based on the projected travel time t\ and mapping (350) the input sample 

15 to the pseudo-offset. In alternate embodiments, the projected source location S' and/or 
the projected receiver location R' are unmoved, or essentially equivalent, to the input 
source location S and the input receiver location R. In even further alternate 
embodiments, the projected source location S' and the projected receiver location R' are 
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projected along the raypath of the input source location S to the scatter point D to the 
input receiver location R. In a further embodiment, the projected source location S' is 
controlled by the ray parameter p? which is solved from equation, 

j 1 ;^ - = h? , (22) 

where the C-wave interval velocity vc(r) is referenced at a datum (401) and d s is the 
distance from the source location S to the datum (401). In an even further embodiment, 
the projected receiver location R' is controlled by the ray parameter ps which is solved 
from equation (21) where d r is the distance from the receiver location R to the datum 
(401). In yet a further embodiment, with p? and p$, the input travel time are computed as 
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(23) 



In still a further embodiment, the projected travel time from S' to D and from D to R' is 
determined by 

1 ^ ro(r) ^ 

trt\ tr<\ «^ 

" I , i+ , , f 1 '*? . (24) 



The input sample point is mapped from / to V . In an even further embodiment, the 
15 pseudo ray parameter p is determined by 

**(/>) = *'. (25) 
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In yet another example embodiment, the pseudo-offset is computed with equation 

(19). 

Those of ordinary skill in the art will recognize that the amplitude term in 
equation (1) contains both the spherical divergence recovery and the obliquity factor 
(Yilmaz, 1987). The filter term in equation (1) represents several functions in alternate 
embodiments. The first function is that of a band-pass filter. The second function is that 
it performs statics correction. The third function is that it operates as a derivative filter, 
which performs the first-order derivative filtering for 3-D data and the half-order 
derivative filtering for 2-D data (Yilmaz, 1987). The last function is that it operates as an 
antialiasing filter. The antialiasing filter is designed for each input/output trace pair. To 
apply the antialiasing for each input/output pair, even the 3-point triangle filter (Lumley 
et al., 1994) will increase the migration runtime significantly. In alternate embodiments, 
to cut the cost and still reject the aliasing noise, an amplitude weighting technique is 
developed. The weighting coefficients are defined as follows, 

4rAr 

where A* is the sample interval, r is the migration radius which is the distance from the 
conversion point (C-wave) or midpoint (P-wave) to the scatter point, Ar is the sample 
interval of the migration radius, and V C i is the C-wave short-spread moveout velocity. 

In another example embodiment of the present invention, a method of migrating 
seismic data is developed for vertical transversely isotropic (VTI) media. The exact travel 
times are implemented for C-wave Kirchhoff prestack time migration: 
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t = t ? + t s . (27) 
For v(z) vertical transverse isotropy (VTI) media, we have the following equations: 



, _ At m v P0i/ v ?i(Pp> v ?Oi> v soi>ej A ) 

P ~^ I 2 2 



(28) 



_ ^ ^sof v so</ v s«(Ps> v Po/» v soi>*i At ) 

/ y~i — » v-v) 

i=l - PSySi(P S ^0i' v S0i^i '3 ) 



(30) 



h / > S V SiO > S» V P0i» V S0i» lg « ) v so,: A (so/ ,, n 

s / — > pi; 

M V 1 -/*^0 , S> v PO/» v S<M'*f A ) 

where A P and h s are the horizontal distance from the source / receiver to the scatter point, 
pe and p s are the P- and S-wave ray parameters, At POi and A? S0 ; are the vertical-travel- 
times in the z'-th layer, vp 0; - and vsot are the vertical velocities of the i-th layer, and vp,- and 

10 v Sl are the group velocities. In a further embodiment, based on the work of Thomsen 
(1986), we derive the group velocities as a function of p P I p s , v POi , v S0! -, £t and S t 
numerically. The ray parameters in equations (28)-(29) are not time-variant, and h P and 
h s are constant for a given scatter point, so we numerically solve equations (30) and (31) 
for/? P and^s- To compute travel times t P and t s , we still need to know Atm, A/soi, v P0/ , 

15 v S o„ Si and S t . Parameters A? S o/ and v SOl are not be measured or estimated directly from P- 
wave and C-wave data. In a further embodiment, we replace them by the C-wave 
vertical-travel-time 
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At 



CO/ 




(32) 



and the interval vertical-travel-time ratio, that is the same as the vertical-velocity-ratio for 
the i-th layer, 



5 In a further embodiment, with equations (32) and (33), At m is also replaced by At C oi 

and yoi> Velocity v P0 / is determined from P-wave short-spread moveout velocity v P2 ^ but 
the P-wave velocity is indexed at P-wave times and is mapped to C-wave times using yoi. 
In a further embodiment, parameter y 0 / is determined by correlating corresponding P- 
wave and C-wave events from stacked sections or near-offsets of prestack gathers and it 
10 is unavoidable with errors. P-wave velocity vp 2 ; often has its own error. If the P-wave 
velocity error is propagated through y 0i - to C-wave migration, it is more difficult to take it 
out without joint processing of both P- and C-wave data. Alternatively, a C-wave 
velocity is used to control the C-wave migration. In still a further embodiment, using the 
definition of the C-wave interval velocity of Thomsen (1999), 



%= A 'so// A '] 



P0/ ~ v P0z7 V S0* * 



(33) 



15 




(34) 



where 



(35) 



4zi =vso;a + 2^.), 



(36) 



= roM - 4 ) , 



(37) 
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we replace v POl - and v S oi with v cu K>u % and 4 Thus, the total diffraction travel time is 
accurately computed with four parameters v a , jt i9 % and % at a given C-wave vertical- 
travel-time fco- 

In various embodiments, we do not have y eff (Thomsen, 1999) which requires P- 
5 wave velocities. The v(z) effect in travel times is reflected in the interval velocity va and 
the interval The VTI effect is reflected in anisotropic parameters s t and S t directly. 

The exact travel time computation does not increase the migration runtime much 
because for prestack time migration, a small travel timetable is built at each scatter point 
and it does not change with the input traces. 
10 In an even further embodiment, the following equations are used to compute the 

pseudo-offset for the C-wave migration of VTI v(z) media 



* = 'p(/>p) + 'sO>s)> (38) 

t*(p) = t, (39) 

t*(p) = h(p) + t s (p) 9 (40) 

x = h v (p) + h s (p), (41) 



Those of ordinary skill in the art will recognize that the pseudo travel time t* is no 
longer hyperbolic. Therefore, velocity analysis based on the pseudo travel time equation 
becomes complicated. 

In an even further embodiment, a method of performing velocity analysis is 
20 provided. As illustrated in Figure 5, the method comprises computing (510) a moveout 
travel time based on an initial model, computing (520) an intermediate travel time, 
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mapping (530) a gather from the moveout travel time to the intermediate travel time, and 
scanning (540) or velocities based on the intermediate travel time. 

In even further alternate embodiments, the gather for velocity analysis is a 
pseudo-offset gather, a migration image gather, a CDP gather, or a CCP gather or any 
other seismic gather that will occur to those of ordinary skill in the art. 

In further alternate embodiments, computing (510) a moveout travel time is based 
on a simple straight ray, a complicated ray tracing algorithm for VTI v(z) media, or any 
other travel time algorithm that will occur to those of ordinary skill in the art. 

In an even further embodiment, computing (520) an intermediate travel time is 
based on an equation which has a scannable parameter for the final velocity. In further 
alternate embodiments, the scannable parameter is a P-wave velocity, an S-wave velocity, 
a C-wave velocity for converted waves or any other scannable parameter that will occur 
to those of ordinary skill in the art. 

In one example embodiment of the present invention, the method for C-wave 
velocity analysis of v(z) VTI media comprises computing the moveout travel time based 
on equation (40) using initial model parameters and computing the intermediate travel 
time based on following equation: 

' 2 °'^('- , 1/ ). A-<B£±. (42) 

where # is the average vertical travel time ratio and V a is a scannable velocity 
parameter. The final C-wave interval velocity v a is determined from V C2 based on the 
following Dix formula: 
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,2 = f CQi V C2 ( f COi ) ~ 1 COM V C2 Q COM ) 
*co; "'COM 



(43) 



In various example embodiments discussed above, travel times, pseudo travel 
times, moveout travel times, projected travel times, have been computed using various 
algorithms. It will occur to those of ordinary skill in the art that travel times used in the 
5 present invention may be determined using any of the various travel time algorithms that 
will occur to those of ordinary skill in the art, 

In a further embodiment, a system is provided for migrating an input seismic data 
point having an input source location and an input receiver location, a scatter point, and 
an image location associated therewith. The system comprises means for determining a 

10 pseudo-offset. In alternate embodiments, means for determining a pseudo-offset 

comprise a computer, a workstation, or software running on any computer, or any other 
means that will occur to those of ordinary skill in the art. The method further comprises 
means for mapping the seismic data point to the image location based at least in part on 
the pseudo-offset. In alternate embodiments, means for mapping comprise a computer, a 

15 workstation, or software running on any computer, or any other means that will occur to 
those of ordinary skill in the art. 

In an even further embodiment, the system further comprises means for 
determining a pseudo-source and a pseudo-receiver. In alternate embodiments, means for 
determining a pseudo-source and a pseudo-receiver comprise a computer, a workstation, 

20 or software running on any computer, or any other means that will occur to those of 
ordinary skill in the art. 
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In an even further embodiment, the system further comprises means for 
determining a pseudo ray parameter. In alternate embodiments, means for determining a 
pseudo ray parameter comprise a computer, a workstation, or software running on any 
computer, or any other means that will occur to those of ordinary skill in the art. 
5 In still a further embodiment, said means for determining a pseudo ray parameter 

further comprises means for determining a pseudo ray parameter dependant at least in 
part upon the ray parameter down from the input source location to the scatter point and 
the ray parameter up from the scatter point to the input receiver location. In alternate 
embodiments, means for determining a pseudo ray parameter comprise a computer, a 

10 workstation, or software running on any computer, or any other means that will occur to 
those of ordinary skill in the art. 

In yet a further embodiment, said means for determining a pseudo-offset further 
comprises means for determining an input travel time associated with the input source 
location, the input receiver location, and the scatter point. In alternate embodiments, 

15 means for determining an input travel time comprise a computer, a workstation, or 
software running on any computer, or any other means that will occur to those of 
ordinary skill in the art. 

In a further embodiment, means for determining a pseudo-offset further comprises 
means for determining a pseudo travel time from the pseudo-source to the pseudo- 

20 receiver via the scatter point. In alternate embodiments, means for determining a pseudo 
travel time comprise a computer, a workstation, or software running on any computer, or 
any other means that will occur to those of ordinary skill in the art. 
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In a further embodiment, the system further comprises a means for migration 
stacking. In alternate embodiments, means for migration stacking comprise a computer, 
a workstation, or software running on any computer, or any other means that will occur to 
those of ordinary skill in the art. 

In a further embodiment, said means for migration stacking further comprises 
means for determining a moveout travel time. In alternate embodiments, means for 
determining a moveout travel time comprise a computer, a workstation, or software 
running on any computer, or any other means that will occur to those of ordinary skill in 
the art. 

In still a further embodiment, the system further comprises means for summing 
the data point along the moveout travel time. In alternate embodiments, means for 
summing comprise a computer, a workstation, or software running on any computer, or 
any other means that will occur to those of ordinary skill in the art. 

In a further embodiment a system is provided for migration of a seismic data point 
having an input source location, an input receiver location, and a scatter point associated 
therewith. The system comprises means for determining a projected source location. In 
alternate embodiments, means for determining a projected source location comprise a 
computer, a workstation, or software running on any computer, or any other means that 
will occur to those of ordinary skill in the art. The system further comprises means for 
determining a projected receiver location. In alternate embodiments, means for 
determining a projected receiver location comprise a computer, a workstation, or 
software running on any computer, or any other means that will occur to those of 
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ordinary skill in the art. The system further comprises means for mapping the seismic 
data point from an input travel time to a projected travel time. In alternate embodiments, 
means for mapping comprise a computer, a workstation, or software running on any 
computer, or any other means that will occur to those of ordinary skill in the art. The 
system further comprises means for determining a pseudo-offset based on the projected 
travel time. In alternate embodiments, means for determining a pseudo-offset comprise a 
computer, a workstation, or software running on any computer, or any other means that 
will occur to those of ordinary skill in the art. The system further comprises means for 
mapping the data point to the pseudo-offset. In alternate embodiments, means for 
mapping comprise a computer, a workstation, or software running on any computer, or 
any other means that will occur to those of ordinary skill in the art. 

In an even further embodiment, the system further comprises means for 
determining a raypath from said input source location to said input receiver location. In 
alternate embodiments, means for determining a raypath comprise a computer, a 
workstation, or software running on any computer, or any other means that will occur to 
those of ordinary skill in the art. 

In still a further embodiment, said means for determining a projected receiver 
location further comprises means for determining a projected receiver location along the 
determined raypath. In alternate embodiments, means for determining a projected 
receiver location along the determined raypath comprise a computer, a workstation, or 
software running on any computer, or any other means that will occur to those of 
ordinary skill in the art. 
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In a further embodiment, said means for determining a projected source location 
further comprises means for determining a projected source location along the 
determined raypath. In alternate embodiments, means for determining a projected source 
location along the determined raypath comprise a computer, a workstation, or software 
5 running on any computer, or any other means that will occur to those of ordinary skill in 
the art. 

In yet a further embodiment, the system further comprises means for determining 
a pseudo-source and a pseudo-receiver. In alternate embodiments, means for determining 
a pseudo-source and a pseudo-receiver comprise a computer, a workstation, or software 
10 running on any computer, or any other means that will occur to those of ordinary skill in 
the art. 

In still a further embodiment, the system further comprises means for determining 
a pseudo ray parameter. In alternate embodiments, means for determining a pseudo ray 
parameter comprise a computer, a workstation, or software running on any computer, or 

15 any other means that will occur to those of ordinary skill in the art. 

In an even further embodiment, said means for determining a pseudo-offset 
further comprises means for determining a projected travel time associated with the 
projected source location, the projected receiver location, and the scatter point. In 
alternate embodiments, means for determining a projected travel time comprise a 

20 computer, a workstation, or software running on any computer, or any other means that 
will occur to those of ordinary skill in the art. 
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In an even further embodiment, means for determining a pseudo-offset further 
comprises means for determining a pseudo travel time from a pseudo-source to a pseudo- 
receiver via the scatter point. In alternate embodiments, means for determining a pseudo 
travel time from a pseudo-source to a pseudo-receiver via the scatter point comprise a 
5 computer, a workstation, or software running on any computer, or any other means that 
will occur to those of ordinary skill in the art. 

In an even further embodiment, the system further comprises means for migration 
stacking. In alternate embodiments, means for migration stacking comprise a computer, 
a workstation, or software running on any computer, or any other means that will occur to 
10 those of ordinary skill in the art. 

In an even further embodiment, said means for migration stacking further 
comprises means for determining a moveout travel time. In alternate embodiments, 
means for determining a moveout travel time comprise a computer, a workstation, or 
software running on any computer, or any other means that will occur to those of 
15 ordinary skill in the art. 

In yet a further embodiment, the system further comprises means for summing the 
data point along the moveout travel time. In alternate embodiments, means for summing 
comprise a computer, a workstation, or software running on any computer, or any other 
means that will occur to those of ordinary skill in the art. 
20 In still a further embodiment, a system for performing velocity analysis on a 

seismic gather is provided. The system comprises means for computing a moveout travel 
time based on an initial model. In alternate embodiments, means for computing a 
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moveout travel time based on an initial model comprise a computer, a workstation, or 
software running on any computer, or any other means that will occur to those of 
ordinary skill in the art. The system further comprises means for computing an 
intermediate travel time. In alternate embodiments, means for computing an intermediate 
travel time comprise a computer, a workstation, or software running on any computer, or 
any other means that will occur to those of ordinary skill in the art. The system further 
comprises means for mapping the gather from the moveout travel time to the intermediate 
travel time. In alternate embodiments, means for mapping comprise a computer, a 
workstation, or software running on any computer, or any other means that will occur to 
those of ordinary skill in the art. The system further comprises means for scanning for 
velocities based on the intermediate travel time. In alternate embodiments, means for 
scanning comprise a computer, a workstation, or software running on any computer, or 
any other means that will occur to those of ordinary skill in the art. 

The example method embodiments are described for the purpose of teaching the 
broad aspect of the invention only. Systems, apparatuses, and the seismic product of the 
methods are also described and claimed as the present invention. For further example 
embodiments of the present invention, readers are directed to exhibits A and B attached, 
incorporated, and made a part of this patent application in compliance with 37 CFR 1.91. 
Exhibit A is an unpublished abstract describing further example embodiments of the 
present invention. Exhibit B is also an unpublished abstract describing further example 
embodiments of the present invention. 
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